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Abstract 

We provide analytic solutions of the nonlinear differential equation system describ- 
ing the particle paths below small-amplitude periodic gravity waves travelling on a 
constant vorticity current. We show that these paths are not closed curves. Some 
solutions can be expressed in terms of Jacobi elliptic functions, others in terms of hy- 
perelliptic functions. We obtain new kinds of particle paths. We make some remarks 
on the stagnation points which could appear in the fluid due to the vorticity. 

1 Introduction 

The present work is confined to two-dimensional water waves with constant 
vorticity. To gain insight into the motion in such waves we determine analyti- 
cally the trajectories beneath small-amplitude water waves in constant vorticity 
flows. In order to get this Lagrangian feature of the flow, that is, the evolution 
of individual water particles, we will firstly find a solution of the Eulerian sys- 
tem of equations within the framework of small amplitude waves and then we 
will calculate the solutions of the nonlinear differential equations system which 
describes the particle motion. 

Many of the theoretical results concerning waves on water make the initial 
assumption of irrotational flow. There are circumstances in which this is well 
justified but there are cases where it is inappropriate. Waves with vorticity 
are commonly seen in nature, for example, in shear currents. Tidal flow is a 
well-known example when constant vorticity flow is an appropriate model (see 
Da Silva T. A. and Peregrine D. H. [16]). For a discussion of the physical 
relevance of flows with constant vorticity see also [S]. In 1802, Gerstner [5D] 
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constructed an explicit example of a periodic travelling wave in water of infinite 
depth with a specific non-constant vorticitj0. The fact that this flow is very 
special is confirmed also by the fact that this is the only steady flow satisfying the 
constraint of constant pressure along the streamlines cf. [S^. Gerstner's wave 
is a two-dimensional wave which adopts the Lagrangian viewpoint, describing 
the evolution of individual water particles. Its surface profile is symmetric [3]. 
Beneath Gerstner's wave it is possible to have a motion of the fluid where all 
particles describe circles with a depth-dependent radius [3], 

In 1934, Dubreil-Jacotin [TB] considered the problem of the existance of 
steady periodic water waves with general vorticity. She proved the existence of 
large classes of small-amplitude water waves with vorticity. For large-amplitude 
water waves with vorticity, in 2004, Constantin and Strauss [l^ proved that, for 
an arbitrary vorticity distribution and for a given c > and relative mass flux pq, 
there is a global continuum of steady periodic waves travelling at speed c in water 
of finite depth and such that the horizontal component of the velocity u < c 
throughout the fluid. The continuum contains waves with u arbitrarily close 
to the wave speed c. The existence of global continua of smooth solutions for 
the related problem of periodic waves of infinite depth was proved by Hur [25] . 
The vorticity does not destroy the symmetry. The construction in [T3] assumes 
that the wave profiles are symmetric. Constantin and Escher [5], [5] proved that 
the symmetry of the wave profile is not a hypothesis but rather a conclusion 
when the wave profile is monotone between crest and trough and the vorticity 
is positive and non-increasing with greater depth. Free from restrictions on 
the vorticity but requiring a quite precise knowledge of all streamlines in the 
fluid, Hur proved in |26j that if the wave profile is monotone near the trough and 
every streamline has a single minimum per wavelength located below the trough, 
then the steady periodic water waves of finite depth are symmetric. Constantin, 
Ehrnstom and Wahlen ffl showed that for an arbitrary vorticity distribution, a 
steady periodic water wave with a profile that is monotone between crests and 
troughs has to be symmetric. 

Another remarkable feature of the rotational steady waves is that they could 
contain stagnation points. If (u, v) denotes the velocity field and c the constant 
speed of the wave, then a point where u — c and w = is called stagnation point. 
There are interesting problems related to the so-called extreme waves: these are 
waves with the stagnation points at their crests. Varvaruca [4^ proved for a 
certain class of vorticity functions, the existence of extreme waves and showed 
that at such a stagnation point the profile of the wave has either a corner of 
120" or a horizontal tangent. For a recent survey of different aspects of the 
theory of steady water waves with vorticity see [39] . 

This paper has interest in finding information about the flow below water 
waves with constant vorticity, more precisely, we will investigate how the pres- 
ence of vorticity influences the particle paths. Throughout the hydrodynamics 
literature, it has been quite common to assume that beneath an irrotational 

^ This solution was independently re-discovered later by Rankine |37) . Modern detailed 
descriptions of this wave are given in the recent papers |3j and |24| . 
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periodic two-dimensional travelling water wave, the particles trace closed, cir- 
cular or elliptic, orbits, (see for example [34], [35], [17], [32]). While in this 
first approximation all particle paths appear to be closed, Constantin and Vil- 
lari showed in [12) . using phase-plane considerations for the nonlinear system 
describing the particle motion, that in linear irrotational periodic gravity water 
waves no particles trajectory is actually closed, unless the free surface is flat. 
Similar results hold for the particle trajectories in irrotational deep-water (see 
Constantin, Ehrnstrom and Villari [B]), and in irrotational shallow water (see 
lonescu-Kruse [57] and [25]) Section 5.1). lonescu-Kruse [37], [2E] obtained the 
exact solutions of the nonlinear differential equation system which describes the 
particle motion in small-amplitude shallow water waves and showed that there 
does not exist a single pattern for all particles: depending on the strength of 
the underling uniform current, some particle trajectories are undulating curves 
to the right, or to the left, others are loops with forward drift, and others are 
not physically acceptable, in the last case it seems necessary to study the full 
nonlinear problem. 

For the full nonlinear problem, Constantin proved in 2]i by analyzing a free 
boundary problem for harmonic functions in a planar domain, that all water 
particles in Stokes waves display a forward drift. For an extension of the in- 
vestigation in [3] to deep-water Stokes waves see Henry [3T]. In a very recent 
paper [TJ], Constantin and Strauss recovered the results in [3] by a simpler 
approach and they also investigated the effect of an underlying current on the 
paths of the particles. While in periodic waves within a period each particle 
experiences a backward- forward motion with a forward drift, Constantin and 
Escher showed in |10j that in a solitary water wave there is no backward mo- 
tion: all particles move in the direction of wave propagation at a positive speed, 
the direction being upwards or downwards if the particle precedes, respectively, 
does not precede the wave crest. 

There have also been some studies of particle paths for rotational waves. 
Within the linear theory, by using phase-plane considerations for the nonlinear 
system describing the particle motion, Ehrnstrom and Villari |19| found that for 
positive constant vorticity, the behavior of the streamlines is the same as for the 
irrotational waves, though the physical particle paths behave differently if the 
size of the vorticity is large enough. For negative vorticity they showed that in a 
frame moving with the wave, the fluid contains a cat's-eye vortex (see [36], Ex. 
2.4). The paper [42] by Wahlen which contains an existence result for small- 
amplitude solutions, based on local bifurcation theory, showed also that the 
predictions for negative vorticity |19j in the linear theory are true. We mention 
that an alternative approach to the existence result in [42 for small-amplitude 
steady waves with constant vorticity was very recently proposed by Constantin 
and Varvaruca [TS]. Beside the phase-plane analysis, the exact solutions of the 
nonlinear system describing the particle motion, allow a better understanding 
of the dynamics. For small-amplitude shallow-water waves with vorticity and 
background flow lonescu-Kruse found in [28l the exact solutions and showed that 
depending on the relation between the initial data and the constant vorticity 
some particles trajectories are undulating curves to the right, or to the left. 
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others are loops with forward drift, or with backward drift, others can follow 
pecuHar shapes (see [28], Fig. 7e). 

Removing the shallow-water restriction, in the present paper we provide ex- 
plicit solutions for the nonlinear system describing the motion of the particles 
beneath small-amplitude gravity waves which propagate on the surface of a con- 
stant vorticity flow. 

In Section 2 we recall the governing equations for gravity water waves. 
In Section 3 we present their nondimensionalisation and scaling. We present two 
different scalings, in one the constant vorticity ujq is scaled whereas in another 
one Wo remains unsealed. We choose x and z the space coordinates, thus, the 
sign of the constant vorticity ujq is opposite to the sign of the constant vorticity 
considered if x and y are chosen the space coordinates. 

In Section 4 we obtain the periodic travelling solutions of the considered lin- 
earized problems (see ([27]), respectively ([lU), and the speed of propagation of 
the linear wave c (see (1^ . respectively The solutions are also written 

in the original physical variables: see ([^ . ([^ . respectively (H5t . pS)) . We 
observe that the speed of the wave and the pressure have different expressions 
in the two linearizations. 

In Section 5 we find the solutions of the nonlinear differential equation systems 
([5T]) and we describe the possible particle trajectories beneath constant vorticity 
water waves. In the study of the system ([5T]) it is interesting to observe that, for 
the first linearization, that is, the one made around still water in which the con- 
stant vorticity wq is scaled, the sign of the wave speed c will infiuence the sign of 
the parameter A which appears in the components u and v of the velocity field. 
Thus, if we consider left-going waves c < 0, we get ^ < and if we consider 
right-going waves c > 0, we obtain A > 0. For the second linearization, that 
is, the one made around a laminar flow characterized by u = ujqz + a, v = 0, 
a being a constant, we obtain that, independent of the sign of ojq, the sign of 
A depends on the sign of c — hotuo — y/ghga, where ho is the finite depth and g 
the constant gravitational acceleration. The expression c — hQOjQ — \/ gh^a could 
be regarded as " the speed" of a wave which can be left-going or right-going. In 
the study of the system ([5T]) a peakon-like trajectory ([7^ comes up (see also 
lonescu-Kruse [3T]). This solution contains the arctanh(-) function, having a 
vertical asymptote in the positive direction (Figure 2). For this solution u — c 
and a stagnation point in the fluid appear only for t ±cxd, where the path 
of the particle has a horizontal tangent. The other solutions of the system (I5ip 
are given by ([ST]) . We show these solutions are not closed curves. Some of these 
solutions can be expressed with the aid of the Jacobi elliptic functions, others 
are expressed with the aid of the hyperelliptic functions. We draw some of the 
curves obtained for different values of the parameters (see Figure 3, Figure 4, 
Figure 5). At the end, we make some remarks on the stagnation points inside 
the fluid. 
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2 The water wave problem 

The two-dimensional gravity waves on constant vorticity water of finite depth 
are described by the following boundary value problem: 



Ut + UUx + VUz 
Vt + UVx + VVz 



-Px 



UJo 



V = r]t + urjx on z = Hq + rj{x, t) 
V — Q on z = 

p = Po on z — ho + ri{x, t) 



(EEs) 

(MC) 
(VE) 

(KBCs) 

(DBC) 



(1) 



where {u{x, z,t),v{x, z,t)) is the velocity field of the water - no motion takes 
place in the y-direction, p{x,z,t) denotes the pressure, g is the constant gravi- 
tational acceleration, po being the constant atmospheric pressure and wq is the 
constant vorticity. The water moves in a domain with a free upper surface at 
z = hQ + r]{x,t), for a constant ho > 0, and a flat bottom at z = 0. We set the 
constant water density p ~ 1- See in the Figure 1 an example of a linear shear 
flow with constant vorticity lj = const := > 0. 




z=0 Bottom 

Figure 1 . Linear shiear fiow. 



3 Non-dimensionalization and scaling 

We non-dimensionalize the set of equations ^ using the undisturbed depth of 
the water ho-, as the vertical scale, a typical wavelength A, as the horizontal 
scale, and a typical amplitude of the surface wave a (for more details see [22]). 
Thus, we define the set of non-dimensional variables 

x I— >■ Ax, z M- hoz, r] I— >■ arj, t —^=t, 

u ^ y/ghou, v^ho^^^v, (2) 



pi-^ Po+ gho{l - z)+ ghop, 
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where, to avoid new notations, we have used the same symbols for the non- 
dimensional variables x, z, rj, t, u, v, p on the right-hand side. The partial 
derivatives Uz and will be then replaced by 

Uz H> — — Uz, H> ho ^„ v^, (3) 
ho V 

and the natural scaling for the vorticity is 

ujo i-> — r — ^^0, (4) 

llo 

where we have used the same symbol for the non-dimensional on the right- 
hand side. 

Therefore, in non-dimensional variables ([2]), (j4|), the water-wave problem ([!]) 
becomes: 

Ut + UU^ + VUz = 
6'^{vt + UV^ +VVz) ^ -Pz 
U^+Vz ^0 

Uz - S'^Voo = Wo (5) 

V — eijit + urjx) on z = 1 + erj{x, t) 

p — erj on z = 1 + eri{x, t) 

V = Q on z = 

where we have introduced the amplitude parameter ^ — and the shallowness 
parameter 6 = 

After the non-dimensionalization of the system ([T]) let us now proceed with 
the scaling transformation. First we observe that, on z = 1 -f ery, both v and p 
are proportional to e. This is consistent with the fact that as e — > we must have 
t; — >■ and p — > 0. We can consider the following scaling of the non-dimensional 
variables 

p 1-^ ep, It i~> eu, V 1-^ ev (6) 

where we avoided again the introduction of a new notation. For this scaling of 
u and V, we also get 

uo f-^ eojo (7) 

The water-wave problem ([1]) writes in non-dimensional scaled variables ([2]), ([4]), 
®, 0, as 

Ut + e{uux + vuz) = -Px 
5'^[vt -f e{uvx + vvz)] = -Pz 
+ = 

Uz - S'^Vx ^ uJo (8) 
V — r/t + turix on z = 1 + eri{x, t) 

p — rj onz=l + eri{x, t) 

V — Q on z = 
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By letting e — >■ 0, 5 being fixed, we obtain a linear approximation of the problem 
([S|), that is, 

Ut+Px =0 

S^vt+p, =0 

Uz - S'^Vx ^ Wo (9) 

V = rjt on z = 1 
p — rj on z = 1 

V = on z = 

This linearization is used in |27j , [281 for irrotational and constant vorticity shal- 
low water waves, in |29j . [30) for capillary-gravity waves and in i31 for constant 
vorticity gravity waves. 

For constant vorticity flows all the streamlines are real-analytic as proved re- 
cently in Thus, for travelling water waves with constant vorticity one can 
get the analytic validity of the linearization ([9]) . In the case of irrotational water 
waves, the rigorous analysis of the validity and relevance of the linearizations 
around some reference states is carried out in [1]. 

Instead of the scaling ^ , we can also consider the following one 

p t-^ ep, u i-^ ujqz + a + eu, v t-^ ev (10) 

with a constant. In this case ujq remains unsealed. Thus, the water-wave 
problem ([T]) writes in non-dimensional scaled variables ([2]), ([4]), (fTO]). as 

ut + e{uux + vuz) + {ujqz + a)ux + ujqv = -Px 
5'^[vt + e{uvx + vvz) + {i^qz + a)vx] = -pz 

Ux+Vz=Q 

Uz - S^Vx = (11) 

V ^ rjt + eurix + f^^omx + (wo + ajrjx on z = 1 + erj{x, t) 

p — r] onz = l + eri{x, t) 

V = on z = 

By letting e —>■ 0, 5 being fixed, we obtain a linear approximation of the problem 
(fTI|) . that is, 

Ut + {ujqz + a)ux + u)qv +Px = 
5'^[vt + {ijJqz + a)vx] +Pz^Q 

Ux+Vz=Q 

Uz - 5\^x = (12) 
1' = ?7t + (^0 + oi)r]x on z = 1 

p = rj on z = 1 

w = on z = 

We observe that the forth equation in the system ([T^ . which represents the 
vorticity equation, becomes in these scaled variables the vorticty equation for 
an irrotational flow. 

The linearization (IT^ is used in for irrotational gravity waves, in [B] for the 
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corresponding deep-water waves, in [22], [23] for capillary-gravity waves and in 
[19] for gravity waves over finite depth with constant vorticity. This linearization 
is around a laminar flow. Such shear flows are characterized by the flat surface, 
z = 1, corresponding to 77 = 0, p = 0, u = and u = ujqz + a. 



4 Solutions of the linearized problems 

Let us solve the linearized systems (|9]) and ([T2|) and compare their solutions. 
From the third equation and the forth equation in we obtain that 

V,, + S\,, = (13) 

Applying the method of separation of variables, we seek the solution of the 
equation (fO)) in the form 

v{x,z,t) ^ F{x,t)G{z,t) (14) 

Substituting pi]) into the equation (IT51) . separating the variables and taking 
into account the expressions of v on the boundaries, that is, the fifth equation 
and the last equation in we find 

v{x, z, t) = . smh{KSz)T]t (15) 

smh(A()) 

where A" > is a constant that might depend on time. Taking into account 
(fTS)) and the fourth equation of the system (|9]), we obtain 

u{x, z, t) = . ■ cos\Y(K6z)T]tx + -I- T(x, t) (16) 

K smh(Ad) 

where J-{x, t) is an arbitrary function. The components u and v of the velocity 
have to fulfill also the third equation in ^ , hence, in view of (USj) and (fTB|) , we 
get 

^ Wr^r X dF{x,t) KS WT^C N 

■ cosh( K 6 z)r]tx^ + — = - . cosh{Kdz)rit (17) 



Ksinh{KS) dx smh{KS) 

The above relation must hold for all values of a; e R, and < z < 1, thus, it 
follows 

= (18) 

ox 

and 

Vt..+K\^0E (19) 
Seeking periodic travelling wave solutions for the equation (jl9p . we take 



K = 27r (20) 



^ For constants K independent on time, we can integrate the equation 1191 1 with respect to 
t and we get 
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and we choose the following solution 

77(x, t) = cos(27r(a; - ct)) (21) 

where c represents the non-dimensional speed of propagation of the linear wave 
and is to be determined. 

From (fTSl) the function J-{x, t) is independent of x, therefore we will denote this 
function by J-{t). Thus, the components of the velocity filed are 

2ttSc 

u(x, z, t) = . , — cos\i(2Ti5z) cos(27r(a; — ct)) + ujqz + F(t) 
smh(2iTd) 

2ttc 

v(x, z, t) = . — — sinh(27r5z) sin(27r(a; — ct)) (22) 
smh(27rc)) 

We return now to the systems ^ in order to find the the expressions of the 
pressure. Taking into account the first two equations in ^ and the expressions 
(22) of the velocity field, we obtain 

2-k6c^ 

p{x, z, t) = — — — cosh(27r5z) cos(27r(a; — ct)) — xF'lt) (23) 

smh(27r(3) 

On the free surface z — 1 the pressure ([25)1 has to fulfill the sixth equation of 
the system Hence, in view of (PTt . we get 

2tt5c^ coth(27r(5) cos(27r(a; - ct)) - xF'{t) = cos(27r(x - ct)) (24) 

The above relation must hold for all values a; G R, therefore, we get 

J-{t) ~ constant :— cq (25) 

and we provide the non-dimensional speed of the linear wave 



2 



tanh(27r(5) 



2tt5 

Summing up, the solution of the linear system ^ is: 
ry(a;, t) — cos(27r(x — ct)) 

Pi^> ^' ^) = sinhf2xa) cosh(27r5z) cos(27r(x - ct)) 

u{x, z, t) = sinhiils) cosh(27rfc) cos(27r(a; - ct)) + ljqz + cq 

^(^^ ^' = sinhr2'7r^) sinh(27r(5z) sin(27r(2; - ct)) 



(26) 



(27) 



The solution of this equation can be written into the form 



r)(x, t) = T{t) [ci cos{Kx) + C2 sin(A'x')] - 



— j R(x) sin{Kx) dx 



cos{Kx) - 



— J R(x) cos{Kx) dx 



sin{Kx) 



Ci, C2 being integration constants. We observe that to the Unearized problem ([9j we can get 
solutions different from the usual one ri(x, t) = cos{K{x — ct)). 



4 Solutions of the linearized problems 



10 



with c given by (1^ . 

Taking into account ([2]), (O, we return to the original physical 

variables. The speed of the wave (|26p and the solution ([?7)) become: 



^— /tanh(fcft.o) / tanh(/c/io) 



77(2;, i) = a cos 



eh{) cos[fc(a; — ct)\ 



(29) 



p{x, z, t)^po+ g{ho - z) + e js^fpls^ cosh(fcz) cos[fc(a; - ct)] 

u{x, z, t) = '^ sinh(khQ) cosh(fcz) cos[fc(a; - ct)] +euoz + e^/ghaCQ 

^(^^ ^' ^) = ^ sinHkho) sinh(fcz) sin[fc(x - ct)] 
where 

^ (30) 

is the wave number. The sign minus in (1281) indicates a left-going wave. 

Let us look now at the linearized system ([T^. From the third equation and 
the forth equation in (|12l) . we obtain again the equation ([T^ . Applying the 
method of separation of variables, we seek the solution of this equation in the 
form (fT4|) . Substituting (fT4|) into the equation (fT3|). separating the variables 
and taking into account the expressions of v on the boundaries, that is, the fifth 
equation and the last equation in (|12p . we find 



v{x, z, t) = . ■ sinh{KSz) [rjt + {ujq + a)r]^] (31) 

smh(iVoj 

where K > is a constant that might depend on time. Taking into account 
(Plj) and the fourth equation of the system (IT^ . we obtain 

u{x, z, t) = . ■ cos\i{K5z) [r^tx + (^^o + oi)t]xx] + ^{x, t) (32) 

A smh(A 0) 

where S'(a;, t) is an arbitrary function. The components u and v of the velocity 
have to fulfill also the third equation in p^ . hence, in view of (|5T|) and ([5^ . 
we get 

■cosh{Kdz)[r]txx + [^^o + a)ria:xx]+ - 



— • u/T^s:\ cosh{KSz) [rjt + (wq + a)r]^] (33) 
smh(Ad) 

The above relation must hold for all values of a; G R, and < z < 1, thus, it 
follows 

^ = (34) 

ox 
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and 

[7]t + [uJa + + K'^ [Vt + i^^a + a)Vx] = (35) 

Seeking periodic travelling wave solutions for the equation psp . we take 

K ^2tt (36) 

and we choose the following solution 

77(x, t) = cos(27r(a; - ct)) (37) 

where c represents the non-dimensional speed of propagation of the linear wave 
and is to be determined. 

From the function S{x, t) is independent of therefore we will denote this 
function by S^(t). Thus, the components (15^ . (|5T|) of the velocity filed are 

u{x, z, t) = ^^"^(^ -^o - a) ^^^^2tt5z) cos{2tt(x ~ ct)) + 'Sit) 
smh(27ro) 

v{x, z, t) = 7/o° ^ sinh(27rfc) sin(27r(a; - ct)) (38) 

sinh(27ro) 

We return now to the systems ((T2|) in order to find the the expressions of the 
pressure. Taking into account the first two equations in ()12p and the expressions 
([38]) of the velocity field, we obtain 

27r(5(c — Wo — a) / N , /„ r \ In I w 

p(x,z,t) = — -7 (c — Wq^: — a) cosh(27rdz) cos(27r(x — ct)) + 

smh(27rd) 

^ ujo{c-ujo- a) ^^^^2Tr6z) cos(2^(x - ct)) ~ xd'it) (39) 
smh(27rd) 

On the free surface z = 1 the pressure ([5^ has to fulfill the sixth equation of 
the system ([12]). Hence, in view of ((37| . we get 



(c — ujQ — a) [2TrS{c — ujq — a) coth(27r5) + ujo] cos(27r(x — ct)) — 

~xS'{t) = cos(27r(a; - ct)) (40) 

The above relation must hold for all values a; € R, therefore, we get 

^{t) = constant := Co (41) 

and the non-dimensional speed of the linear wave c satisfies the relation 

(c - Wo - a) [2nS{c - ujq - a) coth(27r(5) + ujq] ^ 1 (42) 

Solving this equation we find 



—Wo ± + SttS coth(27r(5) 

c = u}o + a-{ . ' .^ (43) 

ind coth(27r(!)) 
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Summing up, the solution of the linear system ([T^ is given by 

ri{x, t) = cos{2it{x — ct)) 

p{x, z, t) = g'^;„^(27r?/ [27r(5(c — ojqz — o) cosh{2iT6z) + 1^0 sinh(27r(52)] cos(27r(a; — ct)) 
u{x, z, t) = ^^^i^i^^^^g "^ cosh(27r<52) cos(27r(x — ct)) + cq 
^' *) = ^''sinh(2'^7r^ sinh(27n52) sm(27r(x - ct)) 



(44) 



with c from (|43p . 

Taking into account ©, Q, pH]) - we return to the original physical variables. 
The speed of the wave (|43|) and the solution p4)) have in physical variables the 
following expressions: 



'gho 



ho -^-0±V^-o' + 4fchocoth(fcho) 
=1^0 + a H 



2khQ coth(kh()) 



houio + oi\/gho + — 



-uiQ tanh(fc/io) lb \J '^'^ tanh'^(fe/io) + 4(;fc tanh(fc/io) 
ri{x, t) = eho cos[fc(a; — ct)] 

p(x, z, t) =po+ g(ho - z) + t ^" ''°'^^(^^.^^^"^ [kho (c - uiqz - yjghoa) cosh{kz) 
+ hdUjQ sinh(fc2)] cos[fc(a; — ct)] 

u(x, z, t) = e Binh{kho)^^"°'"^ cosh(A;2) cos[A;(a; — ct)] + ojqz + a^/ghg + e^/ghoco 

v{x, z, t) = ^ ^-fcoCe^ho^^o^ VgT^oi) gjnh(A,.2) sm[k{x - ct)] 

where 



(45) 



(46) 



(47) 



is the wave number. The solution (HSI) . with a = 0, Cq = was also obtained 
in [H]. 

Comparing with P5)) . we observe that the speed of the wave and 

the pressure have different expressions in the two linearizations. The velocity 
field has in the two linearizations the form: 



u{x, z,t) ^ A cosh(fcz) cos[fc(a; — ct)] + Bz + C 
v{x, z,t) — A sinh(fcz) sin[A:(a; — ct)] 
where, for the linearization 



(48) 



(49) 



khpc 
nh{kho) ' 



B = eujo, C = €y/gh^co 
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and for the linearization 



c - houJQ - a^/gho = 

_ ^ khQ(c-hauiQ-s/ghoa) 
sinh( fc/iQ ) 



—uiQ tanh(A;fco) i ^ cJq ta,nh'^ (kho) + 4gfc tanh(A;/io) 
, B = uJo, C = Q\/gho + e^ghoco 



(50) 



5 Particle trajectories 

Let {x{t), z{t)) be the path of a particle in the fluid domain, with location 
(x(0),z(0)) := {xo,zo) at time t — 0. The motion of the particles below the 
small-amplitude water waves in constant vorticity flows with the velocity field 
is described by the following differential system 



dx 
dt 



u(x, z,t) — A cosh(fcz) cos[fc(a; — ct)] + Bz + C 



^ — v{x, z,t) — A sinh(fcz) sin[fc(a; — ct)] 



(51) 



The values of A, B, C are either (|49|) or ((50|). depending on which lineariza- 
tion we consider. From (j49| . ([50l) . for any ujq ^ 0, we get in the both cases 



A^O 



(52) 



For the first linearization (0), the sign of A depends on the sign of the wave 
speed c. Thus, if we choose in the square root with minus, that is, we 

consider left-going waves, we have A < and if we choose in ^9\ l the square 
root with plus, that is, we consider right- going waves, we get ^ > 0. 

For the second linearization il2\} . the sign of A depends on the sign of 
c — hoUJo ~ \/ghoa. Looking at the expression !i50\) of c — haUJo — y/ghoa, we get 
that independent of the sign of luq, if we choose in i50\) the square root with 
minus then c — hocoQ — ^/ gh^a < 0, thus, A < 0, and if we choose in i50\) the 
square root with plus then c — hoiOQ — \/ gh^a > 0, thus, A > 0. The expression 
c— Hqujo — y/ghoa could be regarded as "the speed" of a wave, which is left-going 
if we take the square root with minus in I150\). and is right-going if we take the 
square root with plus in i50\). 

Indeed, if we choose in (jSOp the square root with minus then, how k and Hq are 
greater then zero, for ujq > 0, the expression — wq tanh(A:/io) — 

\J i-^Q tanh^(fc/io) + 4gfc tanh(/c/io) is evidently smaller then zero, and for 

ujo < 0, —Wo tanh{kho)~ ojq tanh^(A;/io) -|- tanh(A:/io) < is equivalent with 

— Wo tanh(fc/io) < \J i-^o tanh^(fc/io) -I- 4f/fctanh(fc/io)- By raising to the power 
two, the last inequality is equivalent with tanh kho > 0, which is a true inequal- 
ity for k and ho are greater then zero.D 
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To study the exact solution of the system ([?T|) it is more convenient to re- 
write it in the fohowing moving frame 



X — k{x — ct), Z — kz 



This transformation yields 

^ = fcA cosh(Z) cos(X) + BZ + k{C - c) 

^ = kAsmh{Z)sm{X) 
We write the second equation of this system in the form 

dZ 



sinh(Z) 



= k A sin X{t)dt 



Integrating, we get 



If 



log 



tanh 



then 



k A sin X{t)dt 

J kAsinX{t) dt < 

exp k A sin X (t) dt 



Z{t) = 2 arctanh 
Taking into account the formula: 

1 + tanh^ [x 



cosh(2a;) = 



(53) 
(54) 

(55) 

(56) 

(57) 
(58) 

(59) 



1 — tanh^(a;) ' 

and the expression (|58p of Z{t), the first equation of the system (|54l) becomes 



— = kA ^^^ „ cos{X) + 2B arctanh {w) + k(C - c) 
dt \ — 

where we have denoted by 

■w = w{t) := exp kAsinX{t) dt^ 



With ([57|) in view, we have 
From (pTj) we get 



< w < 1 



kAsinX{t) = 



1 dw 
w(Tj~dt 



(60) 

(61) 

(62) 
(63) 
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Differentiating with respect to t this relation, we obtain 



kAcos(X)—— = 
at 

From (1631) we have furthermore 



j2 / , \ 2' 



dw 
Hi 



\dt J 

Thus, taking into account (|64|) . (|65l) . the equation (|60| becomes 

d^w 2w / dw\^ ,9,9 I + 

' - I - k^A^w- - 

1 — 



(64) 



(65) 



dt^ l-w^ \ dt 



dw 



-^k^A'^w^ 
We make the foUowing substitution 

e{w) := k^A^w^ 



2B arctanh (w) + k{C - c) =0 (66) 



dw \ ' 
'dtj 



(67) 



A being different from zero (|5^ . Differentiating with respect to t this relation, 
we get 

^x:: = ^ ^ ^ - ^ (68) 



We replace (|67p . (|68l) into the equation (|66l) and we obtain the equation 
, d^ 2w 



2B arctanh (w) + /c(C - c) $ = 



(69) 



(70) 



' dw ' 1 — w"^ 
A solution of the equation (|69l) is 

e = o 

which, in view of ([ST]) and implies 

sinX(t) = ±l (71) 

Therefore, from ([551) with the condition (|57p . and further from ([55]) . a solution 
of the system i51\) is 



x{t) = ct + consti 

z{t) — J arctanh [exp {—\kAt + const2\)] 



(72) 



consti and const2 are constants determined by the initial conditions (a;(0), z(0)) : 
{xqjZq). This peakon-like solution was also presented in the paper [31]. The 
graph of the parametric curve (1721) is drawn in the Figure 2. 
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Figure 2. Peakon-like trajectory. 

Calculating the derivatives of x{t) and z{t) with respect to t, we get 
x'{t) = c 

2Acxp[— (fcA t+const2)] 



Z'{t) 



l-cxp[-2(kA t+const2)] 

2A cxp[A;A i-f consi2] 
l-cxp[2(feA t+const2)] ' 



, kA t + const2 > 
kA t + const2 < 



(73) 



Hence, for the solution a stagnation point in the fluid, where x'{t) = c, 

z'{t) = 0, appear only for t — > ±oo. We observe that at these points the path of 
the particle has a horizontal tangent. 

The other solutions of the equation satisfy 



dw 



2w 



1 — w'^ 

The homogeneous equation: 



has the solution 



d^ 
dw 



2B arctanh (w) + k{C - c) 



2w 



1 — 



(74) 



(75) 



C(w) = 9{l-w'^) 



(76) 

where is an integration constant. By the method of variation of constants, the 
general solution of the non-homogeneous equation (|74|) is given by 



where 9{w) is a continuous function which satisfies the equation 

1 



de_ 

dw 



2B arctanh (w) + fc(C - c) 



1 — 

The solution of the equation (1751) is 

6{w) = —B arctanh ^(w) — k{C — c) arctanh (w) + /3 



(77) 



(78) 



(79) 



/3 being a constant. Therefore, the solution of the non-homogeneous equation 
((71)) has the expression 

(,{w) = (1 - w^) [/3 - k{C - c) arctanh (w) - B arctanh ^{w)] (80) 



5 Particle trajectories 



17 



Taking into account ([S7|) . we get 



dw 
dt 



(1 - «;2)2 [/3 _ k(C - c) arctanh (w) - B arctanh ^{w)] (81) 

We separate the variables in (|5l]) : 



± ■ 



(1 - w^)^k'^A'^- ^^^_^.yyj, -[13- k{C - c) arctanh (w) - B arctanh '\w)\' 

From (El), (EU), we have 



Z(0 



Thus, 



can be written as 
± 



arctanh {w) 



dZ 



dt (82) 



(83) 



= dt 



tanh^d") 
(l-tanh2(f ))2 



fc(C-e) 



that is, 



± 



dZ 



k^A^ sinh2(Z) - [2/3 - fc(C - c)Z - §Z^] ' 



dt 



By ([55|) we obtain 



X{t) 



1 



dZ(t) 



A:Asinh(Z(i)) dt 
Further, from (j53p . we get another solution of the system i51\) : 

a;(t) ^ ct + i arcsm ^ sinh(Z(t)) 

Z(t) 6em5 the solution of the equation i85\} . 

We observe that the solutions (1871) are not closed curves. 
Indeed, if there exists t2 > ti such that Z{t2) = Z{ti), then, in view of 



(84) 



(85) 



(86) 



(87) 



we 



also have ^(ti) = ±Jk^A^smh^{Z{ti)) - [2/3 - fc(C - c)Z(<i) - §Z^{ti)]' 



f (i2). Thus, 



= ±^k^A^ sinh2(Z(t2)) - [2/3 - k{C - c)Z{t2) - ^Z^t^)] 
although in the moving frame we obtain in this case a closed curve with Z{ti) = 
Z{t2) and X{ti) — X{t2), in the fixed frame we get 2(^2) — z{ti) and ^(^2) — 
x{ti) ~ c{t2 — ti) 7^ 0. If c > 0, the particles which follow these curves will have 
a forward drift, if c < 0, they will have a backward drift.O 
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Let us now investigate more the equation (|85p . Using the formula: sinh (a 

2 



cosh(2z) — 1 ^-j^g equation (j85p can be written in the form 



± = dt 

^ cosh(2Z) - fc^i _ [2/3 - k{C - c)Z - 1^2] 2 

Taking into account the expression of cosh(a;) as Taylor series: 

j,2n 

2! ' 4! ' 6! ' ^ (2^ 

n=0 ^ ' 

we get under the square root in ([M]) the following power series 



X X X X 

cosh(x) = l + - + - + - + ... = ;^-— (89) 



-4/32 + 4fc(C - c)/3Z + [k^A^ + 2/3B - ^^(C - c)^] - Bfc(C - c)Z^ + 



The constant A is different from zero (|52|) . thus, the power series (|90| contains 
for sure powers of Z higher than four. A partial sum of the above series is a 
polynomial of degree higher than four. Thus, considering only a partial sum of 
this series, the solution of the equation (|85p involves a hyperelliptic integral (for 
hyperelliptic integrals see, for example, [5], page 252). Its inversion would lead 
to a hyperelliptic function. 

There are special cases when a hyperelliptic integral can be reduced to an 
elliptic one and thus, its inversion will contain the Jacobi elliptic functions sn, 
cn, sc, etc. If in ()90p we have C = c (by choosing appropriate constants cq, Cq 
in respectively ([50)) ) and we consider powers of Z till six, the solution of 
the equation ([85]) involves the following hyperelliptic integral 



±/ . =t (91) 

2^Z6 + (fe!Ai _ ^) ^4 + (^2^2 + 2/3B) Z2 - 4/32 

We consider the substitution 

Z2 = 1 (92) 
Z 

and thus, the left-hand side in (|9T|) reduces to an elliptic integral of the first 
kind: 

±f , =t (93) 

-2y'-4/32Z3 + (fc2A2 + 2/3B) ^2 + _ ^) z + 

This elliptic integral of the first kind may by reduced to the Legendre normal 
form. 

Case 1 : all the zeroes of the cubic polynomial under the square root in (|93p 
are real and distinct. We denote them by Zi < Z2 < Z3. Because the leading 
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coefficient of tliis cubic polynomial is smaller then zero and its constant term is 
greater then zero, we have either 

< Zi < Z2 < Zs (94) 

or 

Zi < Z2 <0 < Z3 (95) 

Case la: the condition (|94l) is fulfilled. 
Then we introduce the variable (p by (see [38] Ch. VI, §4, page 602) 

Z = Z2 sin^ (f + Z-s cos^ (p > (96) 

and we get 

-4/?2(Z-Zi)(i-Z2)(i-i3) - 

= 4/32 sinVcos^ V{Z3 - ^2)^(^3 - Zi) (l - kf sin^ (fi) > 
dZ = —2 s'mip cos ip{Z3 — Z2)dip 

where the constant < fcj < 1 is given by 

kf := (97) 

Z^-Z, 

Therefore we obtain the Legendre normal form of the integral in (1931) : 

(98) 

the constant factor in front of the integral being equal to 



Ci -.^ ±m^J Z3 - Zi (99) 

The inverse of the integral in (|M)) is sn (the Jacobi elliptic function sine ampli- 
tude, see, for example, [2]) 

sn (Cit;fci) = sinip (100) 
In view of the notations (|92|) and ((96)) we get 

Zit) = ^ — (101) 

^Z2 sn 2 (Ci t; fci) + Z3 cn 2 (Ci fci) 

cn being the Jacobi elliptic function cosine amplitude (see, for example, [2]). 
Taking into account the expressions for the derivatives of sn and cn (see, for 
example, |2]), that is, 

— sn (t; k) — cn (t; k) dn (t] k) 
dt 

— cn (t; k) — — sn (t; k) dn [t; k), 
dt 
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where dn (i; k) := \J\ — k"^ sn ^(t; fc), we obtain 

dZ{t) CiiZ3-Z2)sn (Cit;fci) en (Cit;fci) dn {Cit;ki) 



dt 



Z2 sn 2 (Ci fci) + en 2 (Ci fci) 



(102) 



We introduce pUTj) and in ((ST]) and we get a;(t) and z{t) explicitly. 

We remark that, if ^ - ^ < and k^A^+2l3B > (this can happen, for 
example, for a small enough A and for (B > & /3 > 0) or (S < & /? < 0)), 
the coefficients of the cubic polynomial in (j93p have alternating signs. Thus, by 
Descartes' rule of signs, if all the roots are real, the situation ([M)) occurs. 

Case lb: the condition (jHS)) is fulfilled. 
Zi, Z2, Zj, being the zeroes of the real cubic polynomial under the square root in 
(|93|) . this polynomial has the unique decomposition —4/32(2' — Z\){Z — Z2)(Z — 
Z3). A suitable change of variable transforms the integral (|93|) to the Legendre 
normal form (|98p up to a constant. In general, the elliptic functions can have 
complex arguments (for example, if the constant factor in front of the integral 
(j98p is a complex number, then the obtained sine amplitude function will depend 
on a complex variable) but here we are interested only in the real case. We are 
also looking for a real Z, so, Z introduced by has to be greater then zero. 
With the change of variable , which brings the integral to the Legendre 
normal form (1981) . because now Z2 < 0, we end up with a Z which can be 
positive, negative or zero. Thus, in this case, we get the expression (|10ip of 
Z{t) only if t satisfies 

Z2 sn 2 (Ci t; ki) + Z3 cn 2 (Ci t; ki) > (103) 

that is, by sn 2 (Ci t;ki)+ cn 2 (Ci t;ki) — 1, only if t satisfies 

sn^(Ci^;fci) < . . (104) 
Z3 — Z2 

where < - < 1. 

Z3—Z2 

For a very small positive solution Z3 0, the set of t's which fulfill the above 
inequality (|104p tends to the empty set. Thus, in this case, the hyperelliptic 
integral in (PT|) can not be reduced to an elliptic one and the solution can not 
be expressed with the aid of the Jacobi elliptic functions. The solution will be 
expressed with the aid of a hyperelliptic function obtained by the inversion of 
the integral in (|9T|) . 

Case 2: the cubic polynomial under the square root in (|93|) has only one real 
solution denoted Zq. Because the leading coefficient of this cubic polynomial is 
smaller then zero and its constant term is greater then zero, we have 

Zq>0 (105) 

We denote by p and q the real coefficients such that 

-ip'^Z'^ + ijpA^ + 2/3B) Z'^+i— — \ ZH jg— = -4/32(i-Z„)(i2+pZ+g) (106) 
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We introduce the variable ip by (see [3S] Ch. VI, §4, page 602) 



Z = Zo- ^/Z^+pZo + q tan^ | (107) 
and we get 

-4(3^Z^Zo){Z^+pZ + q) = 

= 4/32 f^Z^+pZo + q) ^-^^ (1 - Sin2 ^) > 

\ / COS r) 



I tan ^ 

dz = -Jzi + pio + g^i^V' 

COS"^ -J 

where the constant < fc| < 1 is given by 

kl -.= 1(1+ ^" + ' I (108) 
V ^Zl+pZ^ + q 

Therefore we obtain the Legendre normal form of the integral in 

1 f dil) 



C2 



t (109) 



'I - kl sin^ -0 

the constant factor in front of the integral being equal to 

C2:=±m{Zl+pZQ + q)-- (110) 
The inverse of the integral in (|109p is 

sn (C2t;k2) = ainijj (111) 
Taking into account (|107p . we get 

m = zo - ^z^,+pz,+q f;';l^\ (112) 

^ 1 + cn (C2 t; k2) 

If t satisfies the following inequality 



'Z^+pZo + q- Zo 
cn {C2 t- k2) > ^ " (113) 

^Z^+pZo + q + Zo 

where -1 < Vfo+P^'>+i- ll < 1 then, Z from (ITT^ is greater than zero and 

y/Z^+pZo+q+Zo 

we obtain by ((92)) the expression of Z{t). We can calculate the time derivative 
of Z{t) and by (|57|) we get x{t) and z{t) explicitly. 
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As in the case lb, for very small positive solution Zo — > 0, the set of fs 
which fulfill the inequality (|113p tends to the empty set. Thus, in this case, 
the hyperelliptic integral in (j91|) can not be reduced to an elliptic one and the 
solution can not be expressed with the aid of the Jacobi elliptic functions. The 
solution will be expressed with the aid of a hyperelliptic function obtained by 
the inversion of the integral in (|9ip . 

We remark that, if - ^ < and k'^A'^ + 2pB < (this can happen, for 
example, for a small enough A and for (B < & ^ > 0) or (B > & ^ < 0)), 

the coefficients of the cubic polynomial in have signs +. Thus, by 

Descartes' rule of signs, only one root is positive and we are in the case lb or 
in the case 2. This positive root it will be close to zero. 

Let us draw below some of the curves obtained for different values of the 
parameters, using Mathematiccd. 

We consider k = 1, ho = I, g — 9.8, e = 0.1, a — 0, (3 — I and ojo = 2 > 0. 
Then, by ([50]) . choosing the square root with sign plus, we get c = 4.07454 > 0, 
A = 0.176526 > 0, B = 2 > 0. We take cq such that C = c. In this case, ah 
the roots of the cubic polynomial under the square root in (|93)) are real and we 
get Z{t) in the form (llOip . The graph of the curve obtained is drawn in Figure 3. 



Figure 3. Undulating curve to the right. 

Using the same values for k, Hq, g, e, a and j3 but taking cjq = 20 > 0, 
we get by (ISHl) . choosing the square root with sign minus, c — 4.29294 > 0, 
A = -1.33654 < 0, S = 20 > 0. We take Cq such that C = c. Then all the roots 
of the cubic polynomial under the square root in ()93|) are real and we get Z{t) 
in the form (llOip . The graph of the curve obtained is depicted in Figure 4. 



Figure 4. Loops with a forward drift. 



For the same k, ho, g, e, a and l3 as above, with cjq = 2 > but choosing the 
square root with sign minus in ([SO)) , we get c = —1.59773 <0, A — —0.306137 < 
0, B = 2 > 0. We take Cq such that C = c. Then all the roots of the cubic 
polynomial under the square root in (1931) are real and we get Z(t) in the form 
(jlOip . The graph of the curve obtained is presented in Figure 5. 



^ In Mathematica the Jacobi elliptic functions are implemented as JacobiSN[n, m := kf]:= 
sn(u;fci), JacobiCN[ti, m := fcj]:=cn(M; fci), JacobiDN[it, m := A;^]:=dn(tt; fci) 
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Figure 5. Loops with a bact(ward drift. 



We choose now wq = —20 < 0, fc, ho, g, e a and /3 having the same val- 
ues as above. We get by (ISO)) , choosing the square root with sign plus, c = 
-4.29294 < 0, v4 = 1.33654 > 0, B = -20 < 0. We take cq such that C = c. 
Then the cubic polynomial under the square root in (|93p has only one real root, 
that is, Zq = 0.000798. The real coefficients p and q from (|106p have the values: 
p = 9.55422, q = 24.8588 and the right hand side in (ITT3| has the value 0.99968. 
We conclude that in this case the solution can not be expressed through Jacobi 
elliptic functions. 



We would like to make some remarks on the stagnation points inside the 
fluid. Calculating the derivatives with respect to t of x{t) and z{t) from ((87|) . 
we get 

1 .inh(Z)^-cosh(Z)(^)^ 



iih(Z)^fe2A2 sinh2(Z)-(4|)^ 



(114) 



1 dZ 
k dt 



X'{t) = 

z'{t) = 

where, taking into account 
(fiZ B 

— = k^A^smh{Z)cos\Y{Z)-[2p~k{C-c)Z-—Z\-k{C-c)-BZ] (115) 

With (|85|) in view, for those Z{t) satisfying the following equation 

B 



kAsmh{Z) 



2P - k{C - c)Z - ^Z^ 



we have 



<PZ 
11^ 



= 0, 



and thus, x'{t), z'{t) from (11141) becomes 

x'{t)^c, z'(t)=0 



(116) 
(117) 
(118) 



Hence, for the solution \81l^ the stagnation points in the fluid are obtained by 
solving the equation U16\) . 

The equation (|116l) can be solved graphically. Depending on the signs and on 
the values of the parameters A, B, C, c and (3, the equation (|116p can have one, 
two, three, four or six solutions. See, for example, in Figure 6 some possibilities 



fcAsinh(Z) 



Which of 



that can occur. With continuous line we have drawn 

these solutions are inside the fluid and their nature can be obtained by a further 
study. 
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